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REGULATION OF THE SPECTRAL PEAK IN GAMMA-RAY BURSTS 
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ABSTRACT 

Observations of gamma-ray bursts (GRBs) indicate that the peak of the burst spectrum forms 
in an opaque region of an ultra-relativistic jet. Recent radiative transfer calculations support this 
picture and show that the spectral peak is inherited from an initially thermal photon distribution 
that has been modified by heating into a non-blackbody spectrum with a high-energy tail. We discuss 
processes that regulate the position of the spectral peak E p and compare the expected range of E p 
with observations. The opaque region of a GRB jet has three radial zones: (1) Planck zone r < Rp 
where a blackbody spectrum is enforced; this zone ends where Thomson optical depth decreases to 
t « 10 5 . (2) Wien zone i?p < r < i?w with Kompaneets parameter y 3> 1 where radiation has a Wicn 
spectrum, and (3) Comptonization zone r > Ryj where radiation spectrum is broadened into a Band 
shape. Besides the initial jet temperature, an important factor regulating E p is internal dissipation 
(of bulk motions and magnetic energy) at large distances from the central engine. Dissipation in the 
Planck zone reduces E p , and dissipation in the Wien zone increases E p . In jets with weak magnetic 



fields, the predicted E p approaches but does not exceed the thermal photon energy near the central 
engine; it is expected to vary around 1 MeV up to a maximum value of about 10 MeV. If the jet carries 
an energetically important magnetic field, E p can be increased by dissipation of magnetic energy in 
the Wien zone. Existing data shows a cutoff in the distribution of E p at about 20 MeV. This suggests 
that Poynting flux dissipation indeed contributes to heating in the Wien zone, at least in the brightest 
bursts with the highest E p . 

Subject headings: plasmas - radiation mechanisms: non-thermal - radiation mechanisms: thermal - 
radiative transfer - scattering - gamma-rays: bursts, theory - relativity 



1. INTRODUCTION 

Observed spectra of gamma-ray bursts (GRBs) peak 
at energy E p that varies around 1 MeV (after correct- 
ing by 1 + z for the cosmological redshift, Kaneko et al. 
2006; Goldstein et al. 2012). The spectrum shape can be 
described by a simple Band function (Band et al. 2009) 
- two power laws that are smoothly connected at E p . 
Bursts of higher luminosity are observed to have higher 
E p . An approximate correlation E p w 0.3 MeV 
was reported (e.g. Ghirlanda et al. 2011), where L 52 
is the burst luminosity (isotropic equivalent) in units of 
10 52 erg s _1 . 

The present paper addresses the origin of the spectral 
peak and processes that regulate its position. 

1.1. Synchrotron model 

A simple phenomenological model posits that we ob- 
serve synchrotron radiation, in analogy with blazar mod- 
els. The model assumes that a nonthermal electron pop- 
ulation is injected in the jet by some dissipative process. 
It gives a spectrum with 



power-law tail, dN^/dE oc iT?/ 2 - 1 at E > E p . 

One possibility for injection of high-energy electrons 
is provided by Fermi acceleration in (mildly relativis- 
tic) internal shocks in the variable jet (Rees & Meszaros 
1994). Recent ab initio (particle- in-cell) simulations of 
collisionless shocks provide conclusive results on the ef- 
ficiency of this process (Sironi & Spitkovsky 2011). The 
post-shock electron distribution is dominated by a quasi- 
Maxwellian component; it peaks at energy 7 pca k?7ieC 2 
that is a significant fraction of proton rest mass m p c 2 , 
7poak = e e (m p /m e ). This gives 
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Here B is the magnetic field measured in the rest frame 
of the jet ("fluid frame"), 7 pC ak is the Lorentz factor at 
which the injected electron distribution peaks (also mea- 
sured in the fluid frame) , and Y is the bulk Lorentz factor 
of the jet itself. If the injection distribution of high- 
energy electrons is a power- law above 7 poa k, dNe/dj oc 



where ri2 is radius in units in 10 12 cm, L52 is the isotropic 
equivalent of the jet power in units of 10 52 erg s _1 , and 
eb is the fraction of jet energy that is carried by the 
magnetic field. If the emission radius was fixed for some 
reason at ~ 10 12 e 2 3 cm then E p would be consistent with 
observations. Significant, perhaps dominant, magnetic 
fields are expected in GRB jets. The field must be ad- 
verted by the conducting plasma from the central en- 
gine, and plausible scenarios (e.g. hyper-accretion disks 
around black holes or proto-neutron stars) invoke strong 
magnetic fields. At radii much larger than the size of the 
central engine the advected field is transverse to the jet 
direction. Simulations of shocks in the presence of trans- 
verse magnetic field with eb > 10~ 4 — 10~ 3 show no 
particle acceleration to high energies 7 3> 7 pca k- Thus, 
synchrotron emission from shocks is not expected to ex- 
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then the synchrotron spectrum has a high-energy tend far above E p 
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The synchrotron model may be viewed more broadly as 
a phenomenological model that does not specify the ori- 
gin of nonthermal electrons. When tested against data, 
the model encounters the following problems. 

(1) Thousands of GRBs have been observed, and most 
of them have E p near 1 MeV (Goldstein et al. 2012). Few 
bursts have E p above 10 MeV and no bursts are known 
with E p > 20 MeV. The synchrotron model does not 
explain the clustering of E p around 1 MeV. The model 
predicts E p oc r7p Cak i? which should give a broad distri- 
bution of E p . There is no reason for this combination of 
B, 7 P oak, and T to be nearly the same in different bursts, 
or even within one burst, as GRBs are strongly variable. 

(2) High-energy electrons quickly cool to 7 < 7 pC ak 
(which makes the process radiatively efficient) and 
should emit radiation at E < E p with photon index 
a = —3/2. A typical low-energy index observed in GRBs 
is a — —1, and many bursts have even harder slopes 
a > (Kaneko et al. 2006). The observed hard slopes 
are in conflict with the synchrotron model. 

(3) The observed spectral peak is sharp. E p may be 
defined as photon energy at which the burst luminosity 
peaks, i.e. where ELe is maximum (Le is the spectral 
luminosity), and the spectrum shape around the maxi- 
mum can be quantified by its width at half maximum, 
E\ < E < E 2 - The observed width \og(Ei/E 2 ) is typ- 
ically 1-1.5 decades in photon energy. The synchrotron 
model predicts a broader peak (see e.g. the predicted 
spectra in Daigne et al. 2011). To make the spectral 
peak as sharp as possible, one has to assume an unreal- 
istic electron distribution that has a step-function cutoff 
at 7 < 7 P eak (Baring & Braby 2004; Burgess et al. 2011). 
It is not reasonable to expect a step-like electron distri- 
bution for a few reasons. First, no realistic acceleration 
process gives electron distribution with a low-energy cut- 
off. Second, a low-energy wing of the distribution must 
be created by fast electron cooling. Third, GRB light 
curves are highly variable, and the expected variations 
in 7 P oak should smear out the cutoff. Note also that syn- 
chrotron spectra with sharp peaks (that could be fitted 
by a Band function) are not observed in any other astro- 
physical objects. A close example is provided by blazar 
spectra (e.g. Ghisellini 2006). Their synchrotron peaks 
(both theoretical and observed) have the half-maximum 
width of several orders of magnitude, much broader than 
in GRBs. 

Similar problems are shared by other versions of opti- 
cally thin emission models, e.g. jitter radiation. Obser- 
vations do not support optically thin models and strongly 
suggest that GRB spectral peaks form in an opaque re- 
gion of the jet. 

1.2. Photo spheric emission 

In the opaque zone of the jet, photons keep interacting 
with the plasma and their spectrum is expected to take 
a well defined shape with a sharp peak. (For example, 
consider the extreme case of a Planck distribution.) The 
radiation is released at the photosphcric radius, and a 
distant observer will see a spectrum with a sharp peak. 

The simplest model of photospheric emission assumes 
a hot, radiation-dominated jet with no baryon loading 
or magnetic field (Paczyhski 1986; Goodman 1986). The 
spectra received by distant observers from such jets are 



exactly Planckian (Beloborodov 2011). The peak energy 
is approximately equal to the average photon energy E, 
which remains constant in the ideal radiation-dominated 
outflow, E = Eq. It may be expressed in terms of the 
jet power Lq and the initial radius ro (comparable to the 
size of the central compact object), 

E p ^E ~ 10 L^r-^ 2 MeV. (3) 

The predicted E p is consistent with observed E p ~ 
0.3 Lg 2 4 MeV if r Q > 10 7 cm and the flow is colli- 
mated within an angle ~ 0.03 ro,7 so that the jet 
power Lq required for apparent luminosity L is reduced, 
Lq ?s 9%L. The Planck spectrum is, however, inconsis- 
tent with data, which rules out the baryon-free, non- 
magnetized jet model. 

GRB jets may be expected to carry baryons and mag- 
netic fields. In this case, the predicted photospheric spec- 
trum is not Planckian. Two types of jets may be consid- 
ered: 

(1) Thermally dominated bary on-loaded jets: at small 
radii the jet is dominated by thermal energy (mainly ra- 
diation). The expanding jet cools adiabatically and its 
energy is converted to the bulk kinetic energy of the bary- 
onic flow (Paczyhski 1990; Shemi & Piran 1990). Any 
subphotospheric heating is expected to change the spec- 
trum emitted at the photosphere (e.g. Rees & Meszaros 
2005; Pe'er et al. 2006). In particular, collisional dissipa- 
tion was shown to peak at optical depths r ~ 10 and its 
detailed calculations yielded spectra consistent with ob- 
servations (Beloborodov 2010; Vurm et al. 2011). Non- 
thermal synchrotron emission was found to significantly 
contribute to photospheric emission at E < E p but never 
dominate the spectral peak. 

(2) Magnetically dominated jets: at small radii the jet 
energy is dominated by the magnetic field. The magnetic 
field gradually dissipates as the jet expands, and the re- 
leased energy converts to the flow bulk kinetic energy. 
The Lorentz factor of the dissipative jet grows with ra- 
dius until more than half of the magnetic energy has been 
dissipated. Drenkhahn & Spruit (2002) studied a specific 
scenario assuming an alternating magnetic field that dis- 
sipates via reconnection. The jet also carries baryons, 
and the ultimate result of dissipation is the conversion of 
magnetic energy to the bulk kinetic energy of baryons. 
This gradual conversion takes several decades in radius; 
it is accomplished at a radius i?diss where dissipation 
peaks and ends. Emitted spectra that may be expected 
in this model have been calculated by Giannios (2008). 

Radiative transfer in a heated subphotospheric region 
has been studied with four different numerical codes 
(Pe'er et al. 2006; Giannios 2008; Beloborodov 2010; 
Vurm et al. 2011), consistently giving a Band-type spec- 
trum. The spectrum broadening from thermal (Bose- 
Einstein) to Band shape occurs at optical depths r 30. 
The resulting E p at r ~ 1 is approximately the same as 
at r - 30 (Beloborodov 2010; Giannios 2012)0 Thus, E p 

1 Radiative transfer simulations show that, without subphotop- 
sheric heating, adiabatic cooling would reduce E v by a factor of 
2t~ 2 / 3 (Beloborodov 2011). Heating offsets this effect and keeps 
E p Ri const at t ^ 30. It does not shift the spectral break to 
higher energies; instead the generated heat at r ^ 30 is spent to 
create a high-energy tail of the spectrum. 
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is mainly regulated by processes at smaller radii where 
t > 30. We focus in this paper on radiative processes 
that occur in this highly opaque zone. 

2. PLANCK AND WIEN ZONES 

2.1. Notation 

We first introduce notation for basic quantities that 
will be used in this paper. Let dL/dQ be the jet power 
per unit solid angle; the corresponding isotropic equiva- 
lent is defined by L = 4ir (dL/dCl). The power is carried 
by photons, protons, neutrons, electrons, pairs, and 
magnetic field. The baryonic component includes pro- 
tons and neutrons; the neutron-to-proton ratio depends 
on the details of the central engine (Beloborodov 2003). 
An important parameter of the jet is its energy per pro- 
ton, 

1=-^-, (4) 

Mc 2 V ' 

where M — An (dM/dCl) is the proton mass outflow rate 
(isotropic equivalent). The proton number density in the 
rest frame of the jet ("fluid frame") is given by 



M 



(5) 



where T(r) is the Lorentz factor of the jet. In addition 
to proton-electron plasma the jet may contain pairs 
of density n±. 

The expansion timescale measured in the fluid frame 
is given by 

T 

tcxp = (6) 

and the characteristic optical depth of the jet at a radius 
r is 

(n + n±) a^r 



r 



(7) 



We will focus on the opaque region t>1. 

The radiation component of the jet power (isotropic 
equivalent) can be expressed as 



L 7 = 47rr 2 -[/ 7 r 2 c, 



(8) 



where is the radiation energy density in the fluid 
frame, and (4/3) U 1 is the radiation enthalpy. The frac- 
tion of the total jet power that is carried by radiation is 
given by 

T A rTT 

(9) 
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As the jet expands, e may evolve as a result of adiabatic 
cooling and dissipative heating. 

2.2. Planck zone 

The central engines of GRBs are hot and fill their jets 
with blackbody radiation of a high temperature. As the 
jet expands, the Planck spectrum is initially enforced 
by huge rates of photon emission and absorption. The 
number density of Planck photons at temperature T is 
given by 

0.2 „ kT 
— 6 3 , 6- -, (10) 



where A = h/m e c is Compton wavelength. The temper- 
ature of blackbody radiation is determined by its energy 
density, aT 4 = U 1 . 

There is a characteristic "Planck radius" Rp inside of 
which radiative processes in the thermal plasma are fast 
enough to enforce a Planck spectrum. Outside Rp pho- 
ton production becomes inefficient and the photon num- 
ber freezes out (until the jet expands to larger radii where 
nonthermal processes may produce additional photons 
via synchrotron radiation). Below we show that the jet 
temperature at Rp is Op re 0.01 — 0.02, almost indepen- 
dent of the GRB parameters. 

Two radiative processes should be considered: 
bremsstrahlung e + p — > e + p + 7 and double Comp- 
ton e + 7 — > e + 7 + 7. Let ns and hp>c be the rates 
of photon production by these two processes [cm -3 s~ 4 ]. 
In a thermal plasma with approximately Planckian radi- 
ation the rates are given by 



(11) 
(12) 



where the numerical factors £ re 0.06 and ip re 0.1 
weakly (logarithmically) depend on O (see Appendix). 
When evaluating the rates of photon production we ne- 
glected the presence of e ± pairs; this assumption will be 
verified below. The total rate of photon production is 

% = "DC + TIB- 

Ratio udc/^b ~ 2(n 7 /n)0 5 / 2 is sensitive to tempera- 
ture and also depends on the photon-to-baryon ratio, 
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Uisng n-f/n ~ 10 5 , one concludes that Hdc > for 
9 <; 0.010 For the estimates of the boundary of the 
Planck zone we will use the approximation n 7 re jidc- 

The balance between emission and absorption of 
Planck photons is maintained as long n-y, 
and radius Rp is defined by the following condition, 

^exp^7 — Tlry. (14) 



Substituting Equations (J6J) and (|T2j). and using Equa- 
tion ||7J) with n± <C n, we can rewrite the condition (| 14[) 

as 

e 2 r = V _1 - (15) 

As radiation is still approximately blackbody at r ~ Rp, 
we can use aT 4 re XJ 1 and substitute U 1 from Equa- 
tion ©. This gives, 



T = 



3 er/n m p & 
4Ta 
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(16) 



From Equations (fT5j) and (fH))) , we find O at the Planck 
radius, 

(17) 



O P -1.3xl0- 2 i?p 1 1 / V/ 6 ?72 1/ 5 6 . 



The scattering optical depth at the Planck radius is rp 



(ipQp) 1 ~ 10 5 . Note that 6p and rp weakly depend 



np 



2 This conclusion may not hold if the jet is polluted with 
clumps or sheets with a low entropy per baryon. In such clumps 
bremsstrahlung could dominate photon production. 
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on the exact value of Rp, which can be determined in 
concrete jet models discussed in Sections 3 and 4. 

Our calculation of 6p neglected e ± pairs. This approx- 
imation is easy to check. At r Si Rp pair annihilation 
balance is maintained with approximately blackbody ra- 
diation field. At <C 1 it gives the following equation 
for positron density n + (e.g. Svensson 1984), 



8C 2 (3) 



n+ (n H 



9 J exp 



(18) 



where £(3) ~ 1.2 is the Riemann zeta function. One can 
see that n + <C n at 6 = Op. 

2.3. Wien zone 

Compton scattering maintains a Bose-Einstein distri- 
bution of photons as long as Kompaneets parameter 



y = 4©r > 1. 



(19) 



Comptonization with y > 1 enforces a common electron- 
photon temperature, T e = T 7 = T . The strong ther- 
mal coupling between photons and electrons may also 
be viewed from the electron point of view. The con- 
dition y ^> 1 is equivalent to tc <C t cxp , where tc = 
3m e c/8aTU-y is the timescale for electron temperature 
relaxation to the radiation temperature. 

The Bose-Einstein distribution is described by the pho- 
ton occupation number 



1 



exp(/i + x) — 1 
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(20) 



where chemical potential fi > describes the deficit of 
photon density compared with the Planck value np. In 
the Planck zone, the distribution takes the Planck form 
(/i< 1) with — np. Outside the Planck radius the 
photon number freezes out and, in a heated jet, n 7 /np < 
1. Thermalized radiation with n 7 <C np is described 
by > 1; then radiation has a Wien spectrum Af ~ 
exp(— /j — x). 

Hereafter we call the region y ^> 1 outside Rp "Wien 
zone" to emphasize the possibility of /i ^> 1, although 
this condition is not required; more exactly, photons are 
required to have a general Bose-Einstein distribution. 
The average photon energy in the fluid frame is between 
2.7kT (Planck) and 3kT (Wien). A reasonable definition 
of the Wien zone r < i?w is given by the condition y > 1 . 

Radiation energy density in the Wien zone is given by 



{J- = 3fcTn 7 , r < Ry/. 



(21) 



As the photon density n 7 may be below n P = aT 4 /2.7kT, 
the radiation density U 1 may be below the blackbody 
value aT 4 . 

Equations (|9]) and (|21[) give the following expression 
for the photon-to-baryon ratio, 



n 7 n nip e 
n 



4m e re' 



(22) 



If radiation were suddenly released to_a distant observer, 
the observed average photon energy E and E p would be 
close to 4rfcT, which gives 



n m p c e 
n-y/n 



(23) 



As known from radiative transfer calculations, the posi- 
tion of the spectral break at E p is weakly changed outside 
i?w- Thus, the observed break must be determined near 
the Wien radius r ~ i?w- Equation ([23]) shows that be- 
sides the parameter n the observed E p is controlled by e 
and the photon-to-baryon ratio in the Wien zone. 

The optical depth at the Wien radius may be expressed 
in terms of E p and T\y = r(i?w) using E p 4rfcT and 

y » i, 

r w m e c 2 

t W ~ s • (24) 

E p 

2.4. Evolution equation for photon-to-baryon ratio 

The transition between the Planck and Wien zones is 
more accurately described by an evolution equation for 
n 7 /n. Using dt — dr/Tc (where t is the proper time 
measured in the fluid frame) the rate of change of n 7 /n 
may be written as 



d_ 

dr 



ft) 



nT c 



(25) 



where n 7 is the net rate of photon production measured 
in the fluid frame, 



n 7 = (n B + n DC ) 1 - 



Up 



(26) 



Here np(O) is the Planck density (Equation [10]), and the 
factor (1 — n 7 /np) takes into account photon absorption; 
the net photon production rate vanishes when n 7 = np 
as expected for blackbody radiation. 

To obtain the equation for function f(r) — n 7 /n first 
note that temperature is related to / by Equation 
which we rewrite as 



e = 



H 



H(r) = 



nnip e(r) 



/ ' " v ' 4m e T(r) ' 
Then from Equation (|25[) we obtain, 



df ( f f 1 ' 2 , , i? 2 



1 - 



AVf 4 
0.2 H 3 



(27) 



(28) 



This equation is easily solved for f(r) for any concrete 
jet model with given T(r), n(r), and heating history e(r). 
The solution also determines O(r) = H/f. 

2.5. Nonthermal processes 

We focus in this paper on the early, opaque stages 
of jet expansion and consider dissipation that produces 
thermal heating of the jet. Nonthermal processes are, 
however, possible even at very high optical depths. Col- 
lisionless shocks or magnetic reconnection could impul- 
sively heat electrons to high energies. Synchrotron emis- 
sion from the accelerated electrons produces additional 
photons that could contribute to n 7 . 

A possible effect of synchrotron emission may be es- 
timated as follows. Suppose a fraction e nt h of the jet 
power is given to accelerated electrons. They immedi- 
ately radiate this power via inverse Compton scattering 
and synchrotron losses. Scattering does not change pho- 
ton number, so only the synchrotron fraction of energy 
losses is relevant; it is given by Ub/U^ — e^/e where 
Ub = B 2 /8ir, €b and e are the fractions of jet energy 



Spectral peak of gamma-ray bursts 



5 



carried by the magnetic field and radiation, respectively. 
The power of synchrotron emission is 



T - €B T 

J-'a — Cnth Li- 

£ 



(29) 



This power peaks at energy E s given by Equation (TTJ), 
assuming 7 pea k is high enough to avoid synchrotron self- 
absorption. The produced photon number, however, 
peaks at the low end of the synchrotron spectrum. The 
photon production may be estimated at0 



h. { Es \ 



1/2 



(30) 



where -E a bs is the energy above which soft photons are 
Comptonized to the Wien peak faster than absorbed. 
Photons may be self-absorbed by the high-energy elec- 
trons or absorbed by the thermal plasma via inverse dou- 
ble Compton effect (Appendix); typically -E a bs ~ 10~ 2 E. 
The synchrotron peak energy E s may be comparable to 
E but it is model dependent. 

Comparing N s with the existing flux of thermal pho- 
tons, one finds 



N~, 



entheB E ( E s \ 1/2 

£ 2 -Babs V ^"abs / 



(31) 



Synchrotron radiation should not affect the photon num- 
ber carried by the jet when this ratio is smaller than 
unity. The highest N s is achieved if E s ~ i? a bs, i.e. 
the accelerated electrons have Lorentz factors just high 
enough to experience synchrotron losses, avoiding self- 
absorption. Even this maximum rate may be insufficient 
to significantly affect the photon number in the jet. 

3. THERMALLY DOMINATED JETS 
3.1. Non-dissipative jet 

Early works on GRBs studied the dynamics of ideal 
(non-dissipative) relativistic hot outflows loaded with 
baryons (e.g. Paczynski 1990; Shemi & Piran 1990). 
The flow acceleration is controlled by parameter n (Equa- 
tion!?]). In a radial flow, the fluid Lorentz factor T grows 
linearly with r until it approaches its maximum value 
r = n at R s ~ nr . 

We consider here jets that remain opaque through R s 
to a larger photosphcric radius R* where r = 1, 



Ler T (l +n±/n) 
Anr rripcr] 3 



(32) 



In the absence of dissipation, the density of relict pairs 
is negligible at R+, n±/n <C 1. 

In a non-dissipative jet, the evolution of radiation is 
fully controlled by adiabatic cooling. Entropy is domi- 
nated by radiation and proportional to the photon num- 
ber; thus, adiabatic cooling conserves photon number. 
The photon-to-baryon ratio n 7 /n remains constant as 

3 The power index of (-Es/Sabs) 1 '' 2 in this equation comes from 
the standard synchrotron spectrum of fast-cooling electrons; it is 
slightly changed if the injected electron distribution is modified by 
e cascade accompanying inverse Compton scattering. 



the jet expands|3 it is set by the initial temperature To 
or the initial radius rn, 



/ & 7 
n 



(33) 



At radii r < R s , where the jet power is dominated by 
radiation, the constancy of photon energy flux L 7 L 

and photon number flux iV 7 (isotropic equivalent) im- 
plies constant energy per photon E(r) — L/N^ — const. 
Adiabatic cooling of photons in the fluid frame is com- 
pensated by the increasing Doppler shift as the jet accel- 
erates, Tr = const. 
At radii R s < r < R+, radiation continues to cool 



adiabatically L T oc n 



1/3 



oc r 



~ 2 / 3 while r = const 



As a result, E is reduced between R s and R* by a factor 
of ~ (i?*/i? s ) -2 / 3 . The resulting observed E p of a non- 
dissipative jet is 



E„ ~ E 



^"^'"rtf MeV. 



(34) 



This standard estimate (e.g. Paczynski 1990) is refined 
by a factor of 2 by accurate radiative transfer calcula- 
tions (Beloborodov 2011). The predicted photosphcric 
emission « (E+/Eq)L is bright and has a high 

E p > 1 MeV only if the jet has n > 600L 



5/32 
52 



3.2. Dissipative jet 

Dissipation can offset adiabatic cooling between R s and 
R+ and keep E close to its initial value Eq. Heating must 
begin somewhere in the Wien zone, at r < Radia- 
tive transfer simulations show that photospheric spec- 
trum suffers an irreversible damage if heating is weak at 
r_< i?w_ an d E at the end of the Wien zone is cooled to 
£\y *C E . Even though strong heating at i?w < r < R+ 
still can increase E to_ ~ Eq , _the photosphcric spectrum 
preserves a break at E\v *C Eq and responds to heating 
by developing a hard slope above Ew . Then the emitted 
spectrum significantly deviates from the Band shape (see 
Fig. 7 in Beloborodov 2010). The early start of heating 
is particularly important for bursts with R+ ~^>_R S where 
adiabatic cooling threatens to greatly reduce Eyj. 

Such deep subphotospheric heating is expected, in par- 
ticular for jets with R+ ^> R s , because they have mod- 
erate Lorentz factors T and thus should start dissipate 
early. For instance, internal shocks develop at radii 
r ~ r 2 A -C i?*, where A is the minimum scale of the 
Lorentz factor fluctuations, possibly comparable to the 
size of the central compact object ~ 10 6 cm. Energy dis- 
sipated in shocks is mostly given to protons, however in 
the region of large optical depth t > 10 this energy is 
quickly passed to electrons via Coulomb collisions (Be- 
loborodov 2010). In effect, shocks produce gradual ther- 
mal electron heating. Compton scattering immediately 
passes the heat received by electrons to radiation, which 
dominates the heat capacity of the jet by a huge factor 
~ riry/n. 

A similar effect is expected in the presence of any other 
mechanism that stirs protons and gives them random 

4 More exactly, the photon-to-baryon ratio is n 7 /(ji + ra n ) where 
n n is the number density of neutrons. We assume for simplicity 
a constant proton fraction Y e = n/(n + n n ) and use n 7 /n as a 
conserved "photon-to-baryon ratio." 
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motions in the fluid frame. It may result from magnetic 
reconnection. Dissipation of subdominant magnetic field 
B 2 /8ir Ss nm p c 2 is sufficient to give thermal energy is Eq 
per photon and keep E from falling far below Eq. 

Coulomb electron heating continues at r > i?w where 
it leads to T e > T 7 and Comptonization of Wien radia- 
tion into a Band spectrum (Beloborodov 2010). Coulomb 
heating is a two-body process and its efficiency decreases 
proportionally to optical depth r as the jet approaches 
photosphere. The resulting photospheric spectrum is 
shaped by heating and transfer effects at optical depths 
r - 10. 

Besides the thermal Coulomb heating, the spectrum is 
also affected by high-energy electrons that are injected 
with energies 7m e c 2 ~ m n c 2 « 140 MeV by n-p col- 
lisions. These electrons produce additional high-energy 
photons via inverse Compton scattering and low-energy 
photons via synchrotron emission. Detailed calculations 
of these processes and their effects on the photospheric 
spectrum are found in Beloborodov (2010) and Vurm et 
al. (2011). 

While Comptonization at r > i?w is controlled by 
the details of collisional heating (and possibly by ad- 
ditional electron heating from collisionless shocks and 
magnetic dissipation), the value of E p remains practi- 
cally unchanged from i?w0 This value is given by Equa- 
tion (|23j) . It depends on the photon-to-baryon ratio n^/n 
that freezes out at the Planck radius Rp as discussed in 
Section 2. For thermally dominated jets with typical 
r] ^ 100 and E p ~ 1 MeV, the Planck radius is reached 
before R s , i.e. photon number freezes out at the radia- 
tion dominated stage when the jet is still accelerating. 

The maximum possible E p corresponds to radiative ef- 
ficiency e = 1 and the minimum photon number n T /n. 
As dissipation can only increase n 7 /n from its central 
value (n 7 /n)o, the maximum E p is achieved if n 7 /n = 
(n 7 /n)rj. It corresponds to 

£ p max « Eq, (35) 

where Eq is given by Equation . 

3.3. Effects of collimation 

Achromatic breaks in GRB afterglow light curves are 
often interpreted as evidence for jet beaming, with a typ- 
ical opening angle of 5-10°. Beaming helps explain the 
extremely high apparent luminosities, up to 10 54 erg s _1 
in some GRBs. Beaming must be achieved through a 
collimation process. It is expected from the pressure 
confinement of the jet by the progenitor star or by a 
non-relativistic dense wind from the outer regions of an 
accretion disk. What effect can beaming have on ob- 
served E p l 

If collimation is not accompanied by significant dissi- 
pation, the expanding jet remains approximately adia- 
batic. It can be described as an ideal relativistic flow 
confined by a wall that determines the cross section of 

5 E p is not exactly equal to E but close to it. For any given 
spectral shape, E p is determined by E. In particular, for a Bose- 
Einstein spectrum, E p pa E. Comptonization at r > Ryj creates 
a high-energy tail above E p that somewhat increases E/E p (not 
counting additional synchrotron photons that dominate at E < 
E p , see Vurm et al. 2011). 



the jet S(r) = So(r/ro)^, where r is the radial distance 
along the jet axis and tq is the size of the central en- 
gine. For instance, a parabolic wall gives S(r) oc r, i.e. 
ip = 1, and uncollimated (radial) expansion is described 
by if) = 2. The opening angle of the jet is determined 
by ip and the radius i? C oii where the wall ends and free 
expansion begins (e.g. where the jet escapes the pro- 
genitor star). Between tq and i? C oii the opening angle 
decreases as 9b ~ (ro/r) 1- '''/ 2 . The jet Lorentz factor 
in the collimation funnel grows as T w (r/r )^/ 2 while 
its temperature decreases as T" 1 (as required by con- 
servation laws, see Section 3.1 in Beloborodov 2003 for 
details). Note that in the funnel T9b » 1; it implies a 
marginal causal contact across the jet. A typical beam- 
ing angle 9b ~ 0.1 at i?diss corresponds to T ~ 10 <C 77 
and temperature ~ O.lGo > 8p, i.e. collimation is 
expected to occur in the Planck zone r < Rp. 

Dissipationlcss collimation conserves entropy, and 
hence does not change photon-to-baryon ratio n 7 /n. 
This implies conservation of total photon number car- 
ried by the jet. It also implies that collimation does not 
change the average photon energy, E. Beaming boosts 
the isotropic equivalent of luminosity L 7 fts L and the 

isotropic equivalent of photon flux N y by the same factor 
~ 9^ 2 = (Rcoii/ro) 2 ^. Their ratio E = L 7 /iV 7 remains 
unchanged from its value at rg, E = Eq. Then photo- 
spheric emission has E — eEo as discussed Section 3.2, so 
a radiatively efficient burst with e ~ 1 basically inherits 
E (and hence E p ) from the central region r ~ tq even 
though L may be strongly increased by beaming. 

Next consider dissipative collimation, for instance colli- 
mation accompanied by shocks (e.g. Lazzati et al. 2009). 
Dissipation generates entropy and hence increases n 7 /n. 
Thus, the total pho_ton number carried by the jet is in- 
creased and hence E (jet energy per photon) is reduced. 
In this case, both E and L are affected by collimation, 
which opens the possibility of a correlated variation of 
E p and L (with 9b being the varying parameter). 

Dissipation slows down the jet at r < R co \i (while in- 
creasing its internal energy). It could also change 9b, i.e. 
dissipation may be associated with changing the posi- 
tion of the "wall." Note that 9b must satisfy T9b 1 
for a causal contact across the jet. This condition is 
marginally satisfied for ideal (non-dissipative) collima- 
tion and easily satisfied for dissipative collimation, as 
dissipation reduces T. Thompson et al. (2007) consid- 
ered the possibility that 9b always takes its maximum 
allowed value 9b ~ and pointed out that this gives 
E oc 9 b ~ 1 and hence E oc L 1 / 2 , consistent with observa- 
tions. The model has, however, two difficulties. First, 
it invokes strong variations in 9b, not only from burst 
to burst but also within a single burst (as an extended 
Ep-L~f correlation is observed in individual GRBs, e.g. 
Ghirlanda et al. 2011). Second, the expected normal- 

1/2 

ization constant K in E p = KL 7 is lower than ob- 
served. The model posits that GRBs of various apparent 
luminosities L ~ Lq9^ 2 have approximately the same 
true power Lq ~ 10 50 erg s _1 , which implies a central 

— 1/2 

temperature kT rs 1 r 06 ' MeV. The brightest bursts 
have E p w 3fcXo, not reduced by dissipative collimation. 

— 1/2 

The expected highest E p w 3r 6 MeV is inconsistent 
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with the observed distribution of E p that extends to 15- 
20 MeV. 

In summary, collimation can dramatically increase the 
apparent luminosity of the burst ._ It leaves E either un- 
changed from its central value E or reduces it below 

Eq. 

3.4. Expected range of E p 

The highest value of E p for a thermally dominated jet 
is determined by 



E p ^E ^ 2.7fcT , aT 4 = U , 



(36) 



where Uq is the energy density at the base of the jet, near 
the central engine^ The energy density is related to the 
jet power by 

L = n r 2 U c(3 , (37) 

where fl is the opening solid angle of the flow near the 
central engine, and (3q — vq/c ~ 1 is the initial fluid 
velocity. This gives 



E n 



-1/4 



4iVo,6 /2 MeV. (38) 



If the jet is driven by an accreting black hole, ro is com- 
parable to a few Schwarzschild radii r g = 2GM/c 2 . The 
maximum possible power is comparable to the accretion 
power, 



GM M a 
ro 



3xl0 53 m 



3r 



erg s 



(39) 



where rh is the accretion rate M acc in units of Mq s^ 1 . 
Realistically, the thermal jet power L is smaller than 
L mlH . The only robust heating mechanism at the base of 
the jet is neutrino-antineutrino annihilation; it supplies 
thermal power (Zalamea & Beloborodov 2011), 



M 

3Mp 



-3/2 



erg s 



(40) 



where x ms is the radius of the marginally stable Kep- 
lerian orbit in units or r 9 ; it is determined by the spin 
parameter of the black hole a and varies between x ms = 3 
for a — and x ms = 1/2 for a = 1. Optimistic assump- 
tions regarding the black hole spin give x ms « 1 (which 
corresponds to a = 0.95). In this case, heating peaks in 
a region of radius ro ~ 3r g » 10 6 (M/3Mq) cm. Note 
that the maximum Lo ~ 10 52 erg s _1 corresponds to a 
minimum ro ~ 10 6 cm. 

In summary, the jet power Lq is sensitive to accretion 
rate as well as the black hole spin. It can vary by sev- 
eral orders of magnitude, which implies huge variations 
in observed luminosity (even when the beaming angle 
Ob remains unchanged). There exists, however, a maxi- 
mum jet power and the corresponding maximum achiev- 
able E p . The estimated maximum L gives _E™ ax =10- 
15 MeV. 

6 A more accurate expression for Uo includes the contribution 
from pairs. The corresponding correction to To (the reduction 
by a factor of [4/11] 1 / 4 ) is compensated by an additional photon 
energy that appears in the expanding jet when the e pairs cool 
down and annihilate. The refo re, the simplified estimate for the 
maximum E v in Equation II36I I is approximately valid. 



4. MAGNETICALLY DOMINATED JETS 

It is possible that the energy output of the central en- 
gine is dominated by a Poynting flux, i.e. carried mainly 
by the magnetic field. Then the total energy per baryon 
(cf. Equation 2]) is given by 



V = Vb + Vth = const. 



(41) 



Near the central engine the contribution of the thermal 
power is small, e = n t h/r/ <C 1; it increases at larger radii 
at the expense of the magnetic part iJb/t]. The magnetic 
field is the main reservoir of energy and a large fraction 
of it must dissipate if the model aims to describe a bright, 
radiatively efficient burst. 

Energy dissipated at radii up to Rp is thermalized 
into Planck radiation with luminosity L~ f {Rp) = r/pMc 2 
(where rjp = 77th [-Rp]) and with the average photon en- 
ergy 

E{R P ) = VP m p c 2 [—) . (42) 

At larger radii the photon-to-baryon ratio freezes out (as 
long as the main photon producer is the thermal plasma, 
i.e. nonthermal processes are less eficient). As dissipa- 
tion continues in the Wien zone, E grows proportionally 
to rjth- The observed E p is associated with E at r ~ i?w> 



E„ 



Vw s/d n Vwm p c 

E(K P ) = — — 

VP { n j/n)R T 



(43) 



where i]w = Vth{Rw)- Equation (j43|) shows that E p 
could be very high if the thermal power is nearly zero at 
r < Rp, i.e. the Poynting flux dissipates in the Wien 
zone and practically does not dissipate in the Planck 
zone. Such a jet would experience "photon starvation" 
- it would have a low photon-to-baryon ratio and the 
dissipated power would be carried by a small number of 
photons with a high energy per photon, hence a high E p . 
The possibility of low n 1 /n suggests that E p in magnet- 
ically dominated jets may exceed the maximum E p of 
thermally dominated jets with the same luminosity and 
baryon loading r\. 

More specific estimates can made as follows. Magneti- 
cally dominated jets are expected to start with a modest 
Lorentz factor at the Alfven radius and then gradually 
accelerate due to dissipation (Drenkhahn & Spruit 2002). 
The dissipated Poynting flux is distributed in compara- 
ble amounts between the thermal luminosity of the jet 
L 1 (dominated by radiation) and the kinetic energy flux 
of baryons, L^ n . This means that T] t h ~ T and the ther- 
mal energy density in the fluid frame is approximately 
equal to the proton rest-mass density, 



(44) 



which gives a simple relation between the jet temperature 
and photon-to-baryon ratio, 



n 



3m P O 



(45) 



The freeze-out of n 7 /n implies approximately constant 
temperature O = kT/mc 2 at r > Rp, i.e. the jet tem- 
perature in the Wien zone remains approximately equal 
to Tp. This fact can be understood in a slightly different 
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way: conservation of photon flux in the Wien zone im- 
plies L 7 oc E « 4F kT; then the scaling L 7 « Tkin T 
requires T — const = Tp. 
From Equation (|f 7j) we have 

— Go^nr* <-» 

Substitution of 9 = 9p to Equation (j45]) gives the 
photon-to-baryon ratio for magnetically dominated jets, 

^ w 6 x 10 4 . (47) 
n 

The observed E p is given by 

E p « | T w 3fcT w « 4r w fcTp a 1.5 MeV, 

(48) 

where Tw = r(i?\y) and Tw — T(Rw) ~ Tp is un- 
changed from the temperature at Planck radius Tp. In 
a GRB dominated by photospheric emission, magnetic 
dissipation is expected to be nearly complete at subpho- 
tospheric radii. Therefore, Tw should be close to its 
asymptotic value, Tw ~ f]. 

A dependence of E p on T was generally expected for 
magnetically dominated jets, even in more approximate 
models. Giannios (2012) estimated E p = AY kT assum- 
ing blackbody radiation T = T e ff as long as y S> 1 
(where aT e 4 ff = C/ 7 ). He found E p oc T 4 / 3 ?? 1 / 3 , which 
gives t? 5 / 3 for r « r\. Comparison of Equation (8) in 
Giannios (2012) with our Equation (|4"8")l shows how the 
existence of the Wien zone with T > T e g affects the 
spectral peak. E p is significantly higher, in particular 
for moderate rj ~ 100, and the scaling with n is linear. 
The predicted E p reaches 15 MeV for I\v = 1000; this is 
comparable to the maximum E p expected for thermally 
dominated jets (Section 3.4). 

E p could be increased above the estimate (|48|) if mag- 
netic dissipation is delayed so that heating is strongly 
suppressed at r < Rp. Then U~ <C nm p c 2 at the Planck 
radius, and the jet is "photon-starved." This regime 
would, however, require even stronger assumptions about 
the central engine, which are already quite demanding for 
any magnetically dominated model that neglects thermal 
photons transported by the jet from the central engine. 

In fact, the model described by Equations (|44|) - (|48l) 
already implies a relatively cold central object. It as- 
sumes that the density of thermal radiation advected 
from the center to Rp is smaller than n m p c 2 , i.e. smaller 
than the contribution from magnetic dissipation. Equiv- 
alently, the initial thermal energy per baryon 77th, o must 
become dynamically unimportant, unable to accelerate 
the jet, before the temperature drops to Tp. Only in 
this case the photon number at Rp is controlled by mag- 
netic dissipation rather than by heat advected from the 
central engine. This condition implies an initial tempera- 
ture To ~ VthfiTp w 577th, o keV. Magnetically dominated 
models assume r] t h,o rj, which leads to a strong upper 
bound on Tq. It is inconsistent, for exampe, with the 
collapsar model that invokes a hot accretion disk with 
strong neutrino emission. A relatively cold central engine 
is expected in the proto-magnetar model (e.g. Metzger 
et al. 2011). 



5. DISCUSSION 

Three characteristic radii are important for the forma- 
tion of a photospheric GRB spectrum: 

(1) Planck radius Rp below which radiation is forced 
to have a Planck spectrum. GRB jets have a well-defined 
temperature Op w 0.01 and Thomson optical depth rp ~ 
10 5 at the Planck radius. The photon number carried 
by a heated (dissipative) jet may change at r < Rp. At 
r > Rp thermal heating is unable to change this number. 

(2) Wien radius Rw below which y ^> 1 and radi- 
ation maintains a Wien (or Bose-Einstein) spectrum. 
Compton scattering enforces thermal coupling between 
the plasma and radiation at r < Rw', however, radiation 
density is below the blackbody density. The Wien zone 
Rp < r < Rw occupies an extended range of optical 
depths, T\y < t < rp, where tw Si 30. The observed E p 
forms at Rw and weakly changed at r > Rw- 

(3) Photospheric radius i?* where Thompson optical 
depth t = 1. Radiation is released around i?*. In dissi- 
pative jets, Comptonization at Rw < r < R* changes the 
released spectrum from Wien to Band shape as a result 
of thermal decoupling of plasma and radiation, T e > T 7 . 

Entropy of GRB jets is strongly dominated by radi- 
ation. Entropy of Planck radiation is proportional to 
photon number. Maintaining a Planck spectrum in a 
non-dissipative jet (where entropy is conserved) does not 
require production of photons, as Planck evolution is 
consistent with constant photon number. Therefore, ra- 
diation in non-dissipative jets keeps a Planck spectrum 
(with a decreasing temperature) inside and outside Rp, 
up to the scattering photosphere where the spectrum is 
modified by transfer effects into a multi-Doppler-shifted 
blackbody (Beloborodov 2010; Pe'er & Ryde 2011). 

Evolution of radiation in dissipative jets is different. 
Dissipation generates entropy, so entropy grows as the 
jet expands. Maintaining a Planck spectrum requires 
a growing photon number flux iV 7 , which is possible as 
long as there are fast reactions producing photons. In a 
jet with thermal plasma, photon production is inefficient 
outside Rp ; the photon number flux freezes out and does 
not keep up anymore with the Planck value. This leads 
to the photon deficit n 7 < np and the Wien spectrum 
between Rp and Rw- 

One may define an effective blackbody temperature of 
photospheric radiation by AttR^ctT* = where L+ — 
T 7 (Tk) is the isotropic equivalent of the photospheric 
luminosity. A popular interpretation of E p for photo- 
spheric emission associates it with the effective black- 
body peak, E p « 3kT+. This interpretation is invalid. 
E p is significantly increased above 3fcTt by dissipation in 
the Wien zone. 

Another simple estimate for E p may be obtained if one 
assumes that radiation is blackbody at r ~ i?w (Gian- 
nios 2012). This assumption also significantly underesti- 
mates E p (Section 4). 

E p of photospheric radiation approximately equals the 
average photon energy at Rw that is given by 

Ep^e(Rw)nm p c 2 [ — ) , (49) 

where e = L~J L is the fraction of the jet energy carried 
by radiation. Equation (|49| assumes that the photon 
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number in the spectral peak is set by the radiative pro- 
cesses in an opaque thermal plasma, which are detailed 
in Section 2. The problem resembles radiation evolution 
in the early universe, although GRB flows are much more 
dissipative. In the universe, the observed photospheric 
radiation (cosmic microwave background) has a nearly 
Planck spectrum and provides stringent upper limits on 
possible subphotospheric dissipation. 

Sections 3 and 4 demonstrate the following difference 
between thermally dominated and magnetically domi- 
nated jets. In the thermally dominated limit, 77th 3> Vb, 
the maximum E p is set by the energy density near the 
central engine (Section 3.4). This follows from the fact 
that the thermal energy per baryon in the expanding 
jet never exceeds its initial (central) value, and that the 
photon number per baryon, n 7 /n, may only be increased 
by dissipation. In the opposite, magnetically dominated 
limit rjB ?7th, E p is mainly controlled by the Lorentz 
factor of the jet (Equation BS)l . In this limit, the ther- 
mal output of the central engine is neglected (i.e. the 
central object is assumed to be cold, see Section 4) and 
heat/photons are gradually generated in the expanding 
jet as a result of magnetic dissipation. 

A plausible scenario invokes comparable contributions 
of the thermal energy and magnetic field to the jet power, 
Vth ~ Vb- Then the maximum E p ~ 10 MeV expected 
for the pure thermal jet can be increased by a factor 
~ 2 by dissipation of additional magnetic energy that 
has been transported by the Poynting flux to r > Rp. 
This picture is consistent with the observed distribution 
of E p that cuts off at about 20 MeV. 

The rate of magnetic dissipation is hard to calculate 
from first principles (e.g. reconnection rate depends 
on the field structure in the jet). As long as photo- 
sphere dominates GRB emission, any energetically domi- 
nant magnetic fields should be dissipated to equipartition 
with thermal energy at radii i?diss < R*- The photo- 
spheric radius steeply decreases with T, i?* oc T -3 , while 
i?diss likely increases with V. Therefore, the condition 
i?diss < R* for a Poynting flux-dominated jet in prac- 
tice means i?diss *C R*, as i?diss ~ R* would require a 



fine-tuning of T. 

This suggests that e# is below unity at r ~ R± in GRB 
jets, regardless of their initial magnetization. Dissipation 
of an energetically dominant Poynting flux is possible but 
would be hidden at high optical depths. It can influence 
E p at r < i?W: but probably not dominate heating at 
r ^3 i?* where the Wien spectrum is Comptonized into 
a Band shape. Heating in the Comptonization zone is 
likely dominated by the collisional mechanism that is ex- 
pected to be strongest at moderately subphotospheric 
radii r ^ i?*. It begins at R n ~ (cr n /<7x)-R*, where 
a n w ctt/20 is the nuclear cross section, and dissipates a 
large fraction of the jet energy. 

A positive correlation between E p and the burst lu- 
minosity is expected in the model. Thermally domi- 
nated jets with fixed tq — const and beaming angle 
8b = const would have E p scaling as L 1 / 4 . The E p -L cor- 
relation steepens in the presence of magnetic dissipation 
if brighter bursts have higher T (see also Giannios 2012). 
Recent analysis (Ghirlanda et al. 2011) suggests that the 
slope of the correlation increases as the observed L in- 
creases. This would be consistent with brightest bursts 
having a bigger contribution from the Poynting flux dis- 
sipation. The correlation slope can also be affected by 
dissipative collimation, which may lead to E p cx L 1 / 2 
(Thompson et al. 2007), although collimation alone can 
hardly explain the entire range of the correlation (Sec- 
tion 3.3). 

The theoretical E p discussed in this paper is deter- 
mined not only by the burst luminosity L but also by 
the photon-to-baryon ratio in the jet. This helps under- 
stand the pattern of E p variations in individual pulses of 
GRB light curves. Although a tracking behavior E p (L) 
is generally expected, the presence of the second param- 
eter can lead to significant deviations. In particular, the 
observed high E p at the beginning of a pulse, before its 
luminosity reaches maximum, may be explained as emis- 
sion with a low photon-to-baryon ratio. As the pulse 
progresses, n 7 /n quickly grows until the normal tracking 
behavior is established. 

This work was supported by NSF grant AST-1008334. 



APPENDIX 

A. PHOTON PRODUCTION RATES 

The rate of photon production by the double Compton effect has been extensively discussed in the literature (e.g. 
Thorne 1981; Lightman 1981; Pozdnyakov et al. 1983; Svensson 1984; Chluba et at. 2007). Scattering of mono- 
energetic photons of energy hi>o and density n 7 on cold electrons of density n produces secondary photons hv with the 
following differential rate, 



dfi 



DC 



din; 



4a 
3tt 



n n 7 <7t c 



(Al) 



where a — e 2 /he = 1/137 is the fine structure constant. Equation (|A1[) is valid when the secondary photon has energy 
hv <C hvQ. Extension to v ~ vq was discussed by Gould (1984); it gives only a small correction to the total photon 
production rate 7t.dc, as 77.dc has a flat distribution ov er In v and most photons are emitted with v -C vq. 

For primary photons hv Q with a given spectrum, rate (|A1|) should be averaged over the spectrum. For a Bose-Einstein 
radiation with temperature T, this gives 



dn 



DC 



din: 



4a 
3^ 



n n 7 <7t c Xq O 2 



xq 



hv 
kT' 



hv 
kf' 



e = 



kT 



(A2) 



where x\ 
contribution. 



10.35 for a Planck spectrum and Xq = 12 for a Wien spectrum; photons with xq ij 1 make the dominant 
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The net rate nnc is obtained by integrating Equation (| A2|) over x from a minimum value x m ; n to x ~ 1. Here x m i n 
is a minimum energy of produced photons that avoid absorption and get a chance to be Comptonized to the Wien 
peak (Thorne 1981; Lightman 1981). It is determined by equating the absorption rate t~j~ = «dc c (where odc is 
the absorption coefficient due to inverse double Compton) and the Comptonization rate (where t\c is the time it 
takes to Comptonize the photon energy by a factor of 2). Absorption coefficient aoc is found from Kirchhoff's law 
for Rayleigh- Jeans radiation a^cQitv 2 kT/c 2 — hv dnoc I 'dv : Comptonization rate at non-relativistic temperatures is 
t^Q « AQnaTC. This gives, 



V2 /8 

= -aC(5)9 



1/2 



0.149 1/2 , 



(A3) 



where A = h/m e c is Compton wavelength. We hav e used the relation X 3 n 7 XQ = (24/7r 2 )£(5) O 3 for blackbody radiation, 
where £(5) w 1.037. Integration of Equation (IA2[) gives 



n DC ^ ipnn 1 aTC<d , ^ = v x o ln x min 

07T 



(A4) 



This equation assumes that the scattering electrons are cold in the sense that 8 < 1. In the case of a hot plasma it 
should be multiplied by a correction factor that has been obtained by Svensson (1984) (see also Chluba et al. 2007). 
The correction factor is given by 



5dc = (1 + 13.916 + 11.059 2 + 19.99 3 



(A5) 



At the boundary of the Planck zone we find 0p « 10 2 (Section 2.2); then ip ps 0.1. 

Photon production by bremsstrahlung is given by (e.g. Illarionov & Sunyaev 1975; Thorne 1981; Pozdnyakov et al. 
1983), 



dfiB 
din a; 



3/2 



<\ n^a-r cG 1 ^ 2 ln — . 

x 



(A6) 



Relativistic corrections (e.g. Svensson 1984) are small for temperatures of interest here. Comparison with Equa- 
tion (|A2I) shows that at relevant temperature 9p 0.01 and photon-to-baryon ratio n 7 /n ~ 1 5 , bremsstrahlung 
emissivity is smaller or comparable to the double Compton emissivity. Integration of Equation (|A6[) over ln x from 
In a; m in to ln x ~ gives 



"DC = £,n 2 cr T c6 1/2 , 
Near the boundary of the Planck zone £ ~ 0.06. 



3/2 



"( lna; min) 2 



(A7) 
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